Dispersion of the odd magnetic resonant mode in near-optimally doped 

Bi 2 Sr2CaCu20 8+( 5 
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We report a neutron scattering study of the spin excitation spectrum in the superconducting 
state of slightly overdoped Bi2Sr2CaCu20s+a system (T c =87 K). We focus on the dispersion of the 
resonance peak in the superconducting state that is due to a S=l collective mode. The measured 
spin excitation spectrum bears a strong similarity to the spectrum of the YBa2Cu306+ x system for 
a similar doping level (i.e. x ~ 0.95 — 1), which consists of intersecting upward- and downward- 
dispersing branches. A close comparison of the threshold of the electron-hole spin flip continuum, 
deduced from angle resolved photo-emission measurements in the same system, indicates that the 
magnetic response in the superconducting state is confined, in both energy and momentum, below 
the gapped Stoner continuum. In contrast to YBa2Cu3C>6+x, the spin excitation spectrum is broader 
than the experimental resolution. In the framework of an itinerant-electron model, we quantitatively 
relate this intrinsic energy width to the superconducting gap distribution observed in scanning 
tunnelling microscopy experiments. Our study further suggests a significant in-plane anisotropy of 
the magnetic response. 



I. INTRODUCTION 

It is now well established that the spin excitation spec- 
trum in the superconducting (SC) state of many high- 
T c superconductors is dominated by an unusual spin 
triplet excitationi^ii^ii^^iiiii 2 .. This excitation is 
referred to as the magnetic resonance peak. It is cen- 
tered at the planar antiferromagnetic (AF) wave vector 
Qaf = (tt/o, 7r/a) (where a is the lattice spacing) and at 
an energy E r that scales with the SC critical temperature, 
T c . Further, the resonance peak intensity vanishes above 
T c and exhibits a temperature dependence that looks like 
an order parameter. A renormalization of its character- 
istic energy with temperature is not observed within the 
experimental error^ 6 -. This behavior has been reported in 
several families of copper oxides with maximum SC crit- 
ical temperatures T™ ax >90 K: in Tl 2 Ba 2 Cu0 6+x with 
uniformly spaced, single Cu0 2 layers^ 2 -, as well as in bi- 
layer systems such as YBa 2 Cu 3 6+x (YBCQjW^AW 
and Bi 2 Sr 2 CaCu 2 8+x (Bi2212)i2»ii. The magnetic res- 
onance peak indicates the existence of a S=l collective 
mode with a peculiar dispersion. Its downward dis- 
persion starting at qAF was first observed in YBCO 13 . 
Complementary measurements provided strong indica- 
tions of a second upward dispersion starting from the 
resonance energy in optimally doped YBC Q 14 ' 15 i 16 . Im- 
portantly, both dispersive branches vanish above T c as 
does the magnetic resonance peak at (\af- In strongly 
underdoped cuprates (T c < 62 K), similar dispersive ex- 
citations have also been report e d 17 i 18 ' 19 i 20 ' 21 , but only 
the downward branch vanishes above T c whereas the 
upward dispersion remains essentially unchanged across 
T^. Therefore, the spin excitation dispersion in YBCO 
exhibits a "hour glass" -like shape centered at the reso- 



nance peak. 

Recently, the debate became focused on the origin of 
the S=l dispersive collective mode. The theoretical de- 
scription of the mode is important, because antiferromag- 
netism is generally believed to play a significant role in 
the SC pairing mechanism in high-T c cuprates^. 

First, based on the spin dynamics data in the stripe 
ordered system La7/gBai/gCu04, it has been proposed 
that E r could be a saddle point in the dispersion, with 
spin excitations propagating along a given in-plane di- 
rection (say a*) below E r and along the perpendicular 
direction (i.e. b*) above E, 23 i 24 . This gives rise to an 
X-like shape in twinned crystals where a* and b* are 
mixed. Saddle points can arise in models with spin and 
charge stripe order at low temperature. The spin excita- 
tion spectrum can then be modelled by a specific bond- 
centered stripe model according to which non-magnetic 
charge stripes separate a set of weakly coupled two-leg 
spin ladders in the copper oxide layers. The low en- 
ergy excitations (below E r ) correspond to collective ex- 
citations that propagate perpendicular to the ladder di- 
rection, whereas the high energy part of the spectrum 
(above E r ) is associated to intra-spin ladder excitations 
propagating parallel to the lines of charges. This pic- 
ture, later on sustained by calculation s 25 ! 26 ' 27 , implies 
a pronounced in-plane anisotropy of the magnetic spec- 
trum. However, this is not consistent with the spin exci- 
tations in SC cuprates, in particular in underdoped and 
optimally doped YBCO 20 ' 21 . Indeed, using detwinned 
YBCO samples^, it has been shown that the spin ex- 
citation spectrum exhibits a 2D geometry both below 20 
and above^i E r , inconsistent with a saddle-point disper- 
sion. However, recent calculations of the spin dynamics 
considering fluctuating stripe segments™ shows that the 
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magnetic spectrum is losing its ID character for short 
charge segments (a few atomic distances) because they 
actually exist in both perpendicular directions. Further, 
it should be emphasized, once again, that the resonance 
peak intensity in all SC cuprates exhibits systematically 
a strong temperature dependence in the SC state 1 . This 
is in a marked contrast to the data in the stripe-ordered 
system, where this anomalous temperature dependence 
is absent. These inconsistencies cast some doubt about 
a similar origin of the resonance peak seen in cuprates 
where superconductivity is well developed 1 - and the one 
reported in the stripe-ordered system 2 ^. As a matter 
of fact, it would be more meaningful to compare the 
spectrum in non-SC L^/gBai/gCuG^ with the magnetic 
spectrum in the normal state of other SC cuprates, as it 
has recently been done in underdoped YBCO^i. 

Second, starting from the metallic side of the phase di- 
agram of high-T c superconductors and reducing the hole 
doping, one can try to understand the S=l collective 
mode within an itinerant-electron model. It has been 
proposed that the resonant magnetic collective mode 
could be described as a spin excrfo7 i 29 i 30 ' 31 i 32 ' 33 i 34 , i.e 
a S=l bound state, pushed below the gapped Stoner 
continuum in SC state by AF interactions. This type 
of excitation exists in the SC state only and vanishes 
when the gap disappears in the normal state. Alterna- 
tively, when starting from the Mott-insulator side of the 
phase diagram and increasing the hole doping, in frame- 
work of localized-spin models the mode can be viewed 
as the remnant of the magnon observed in the insu- 
lating AF stat o 35 i 36 . The collective modes of localized 
spins on Cu sites may survive in the metallic state, but 
are heavily damped by scattering from charge carriers. 
Long-lifetime collective excitations can then be restored 
in the SC state, when scattering processes are eliminated 
below the gapped Stoner continuum. When the mode 
energy is close to the gap, it can be viewed as a spin- 
exciton, as in the itinerant-spin approach 35,36 . Both ap- 
proaches represent two different limits of a dual descrip- 
tion of the magnetism of high-T c superconductors: lo- 
calized and itinerant spins are tightly bound and cannot 
be disentangle d 37 ' 38 . It is worth emphasizing that in a 
dual approach one can schematically ascribe the upper 
dispersion to the localized character of the magnetic re- 
sponse and the lower one to the itinerant one. However, 
future quantitative calculations are necessary to validate 
this picture. 

In all of these models, the change of the band elec- 
tronic excitations upon passing through T c has an im- 
portant feedback on the spin excitation spectrum in the 
SC state. The determination of the electron-hole spin 
hip continuum then requires a good knowledge of the 
fermionic dispersion relations, which is still in a stage of 
rapid development in the YBCO system^ 9 .. In contrast, 
the charge excitations in Bi2Sr2CaCu20s+5 have been 
extensively studied by surface-sensitive techniques. First, 
the Fermi surface and the band dispersions have been de- 
termined by angle resolved photo-emission spectroscopy 



(ARPES) (se o 29 i 40 ' 41 and references therein). Second, 
Scanning Tunnelling Microscopy (STM) dat a 42 ' 43 ' 44 ev- 
idence a local distribution of the superconducting gap 
and then infer the low energy quasi-particle excitations 
and the Fermi surface through Friedel oscillations. Re- 
cently, attempts to compute the spin excitation spec- 
trum in nearly optimally doped Bi2212 starting from 
ARPES measurements 4 ^ 6 - have been performed using 
the itinerant-spin approach (excitonic scenario). 

In this paper, we present a study of both energy and 
momentum dependences of the resonant spin excitations 
in the SC state of a nearly optimally doped Bi2212. Our 
study reveals that the spin excitation spectrum in Bi2212 
bears close similarity to the spin excitation spectrum re- 
ported in YBCO for the same doping level. We fur- 
ther show that the resonant spin collective modes are 
located below the gapped Stoner continuum computed 
from ARPES data. The combination of our inelastic 
neutron scattering (INS) measurements and ARPES and 
STM measurements on the same system provides the op- 
portunity to test theoretical models for the spin dynamics 
in the cuprates. 

II. EXPERIMENTAL 

For the present study, we used a single crystal of 
of slightly overdoped Bi2212 with volume ~1.5 g and 
T c =87 K . The INS measurements were performed on 
the thermal triple-axis spectrometers IN8 at the Insti- 
tute Laue Langevin in Grenoble (France), 2T at the Lab- 
oratoire Leon Brillouin at the reactor Orphcc in Saclay 
(France), and PUMA at the reactor FRM-II in Garch- 
ing (Germany). Measurements were carried out with a 
double-focusing PG(002) monochromator and a PG(002) 
analyzer. The final neutron wave vector was set to 
fc/=4.1 A -1 , yielding an energy resolution ~6 meV. 
A PG filter was inserted into the scattered beam in or- 
der to eliminate higher order contaminations. For high 
energy transfers (hu) >55 meV), the PG filter was re- 
moved and kf set to 5.5 A -1 . To cover the full q- 
dependence of the spin excitations, three different scat- 
tering planes were used for the measurements. The sam- 
ple was successively oriented such that momentum trans- 
fers Q of the form (H,H,L), (K/3,K,L) and (H,0.3L,L) 
were accessible. We use a notation in which Q is in- 
dexed in units of the tetragonal reciprocal lattice vectors 
2n/a = 1.64A- 1 = 2n/b and 2n/c = 0.203A- 1 . 

The magnetic neutron scattering cross section is pro- 
portional to Imx(Q,u>), the imaginary part of the dy- 
namical magnetic susceptibility, weighted by the square 
of the Cu magnetic form factor, F(Q), and the detailed 
balance temperature factor— For a paramagnetic sys- 
tem, it reads: 

4n~ = ^\ F (Q)\ 2 - 1 -^—Im X (Q,u J ) (1) 

dUdhj 2n i_exp(--^^) 

where Q = (H, K,L) is the full wave vector and q = 
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FIG. 1: Difference between constant-Q scans performed 
at 10 K and 100 K. INS measurements were carried out 
on spectrometer 2T. The differential spectra are shown at 
different wave vectors of the form: a) Q=(H,H,L=14), b) 
Q=(K/3,K,L=4.2). In addition to the energy scan data 
(open circles) are also reported the magnetic intensities 
(open squares) and the location of the negative background 
(crosses), deduced from the analysis of constant energy scans 
(Fig. [2]). Those data , obtained on spectrometer IN8, have 
been rescaled. The lines indicate the energy dependence of 
the negative background (see text) and the enhancement of 
the magnetic intensity is described by a set of Gaussian func- 
tions on top of the negative background. 



(H, K) is the planar wave vector in the Cu02 plane, 
ro = 0.54 10~ 12 cm is the neutron magnetic scattering 
length. Like the YBCO system, the Bi2212 system con- 
tains two CuC>2 planes per unit cell. Owing to the in- 
teraction between the CuC>2 planes within a bilayer, spin 
excitations with odd (o) and even (e) symmetry with re- 
spect to the exchange of the layer contribute to the spin 
susceptibility as in YBCO^^i, so that: 

x(Q,w) = sin 2 (7rzL)xo(q,w) +cos 2 (7rzi)x e (q,w) (2) 

z=0.109 is the reduced distance between the Cu02 planes 
of the bilayer. Using different L values of the c* compo- 
nent of the momentum transfer, the mode of each sym- 
metry can be measured. The observation of resonant 
modes with both symmetries has recently been reported 
in Bi2212 4 - in two samples and in particular in the sam- 
ple studied here. In the present study, we focus on the 
odd spin excitations, that can be selectively probed by an 
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FIG. 2: Differences between constant energy scans performed 
at 10 K and 100K: a) scans along the (110) direction, b) scans 
along the (130) direction. All measurements were carried out 
on the spectrometer IN8. The position of the negative back- 
ground used in the analysis of the energy scans, reported in 
Figfj] is also shown at the AF wave vector (red circles). The 
solids lines correspond to the fit of the data by a single or 
double Gaussian functional form, on top of a sloping back- 
ground. 



appropriate choice of the L component: L=14 for the pla- 
nar wave vector (0.5,0.5) and L=4.7 for the planar wave 
vector (0.5,1.5). Note that the square of the Cu magnetic 
form factor is about 1.4 times larger for (0.5,1.5,4.7) than 
for (0.5,0.5,14). This is because the magnetic form factor 
is anisotropic in the cuprates^. 



III. INELASTIC NEUTRON SCATTERING 
MEASUREMENTS 

Throughout this paper, we focus on the enhancement 
of the magnetic response in the SC state, corresponding 
to the resonant spin excitations. To extract this response, 
we subtracted scans at T=10 K and T=100 K (> T c =87 
K). This procedure generates a negative background in 
the differential spectra, owing to the thermal enhance- 
ment of the nuclear background. The main contribu- 
tion to this negative background comes from the ther- 
mal population of phonons given by the detailed balance 
factor. In the energy range of interest for the present 
study, 25-60 meV, the nuclear (phononic) background 
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Refs. 


T C (K) E r (meV) 


ov(meV) 


cr(meV) 


Fong et a*ili 


91 


43 


13 ± 2 


12 ± 2 


Present study 


87 


42 


13 ± 2 


11 ± 2 


He et aJii 


83 


38 


12 ± 2 


10 ± 2 


Capogna et a££ 


70 


34 


8 ± 1 


5 ± 1 



TABLE I: Characteristic energy and energy width o> 
(FWHM: Full Width at Half maximum) of the magnetic res- 
onance peak reported until now in Bi2212 samples, a stands 
for the intrinsic FWHM of the resonance peak, after decon- 
volution by the energy resolution ~ 6 meV. 



continuously decreases. Thus, the energy dependent neg- 
ative background can be quite well approximated by a 
functional form like (a + buj) / (exp(— fc ^ ) — 1), with 
T o =100K. Owing to the weakness of the magnetic signal 
( < 10 % of the nuclear one), one also becomes sensi- 
tive to multiple scattering effects. This is the main rea- 
son why the background in the difference signal remains 
weakly negative even at high energy, where the detailed 
balance factor has a negligible effect. Difference scans 
performed far away from the antiferromagnetic wave vec- 
tor at Q=(0.275,0.275,14) (Fig.[T]a) and Q=(0. 5, 1.8,4.7) 
(Fig-fflb) show the typical energy dependence of the neg- 
ative background. On top of the above described negative 
background, the enhanced magnetic intensity appears in 
energy scans at the AF wave vector at Q=(0.5,0.5,14) 
(Fig. Ha) and Q=(0.5, 1.5,4.7) (Fig.fflb). The resonant 
magnetic signal is located at E r =A2 meV and exhibits a 
Gaussian profile with an energy width of a r ~ 13 meV 
(FWHM). From constant-Q scans (Fig. [I]), it appears 
that the magnetic intensity at all wave vectors is limited 
at low energy by a spin-gap of ~ 32 meV. Below the spin 
gap, the magnetic signal is, at least, less than 1/4 of its 
value at E r =42 meV. This value agrees with reports on 
YBCO for similar doping levels*^^. 

In addition to the constant Q-scans shown in Fig. [TJ 
systematic constant-energy scans have been performed 
in each scattering plane (Fig. to characterize the q- 
dependence of the magnetic intensity. Along the (130) 
direction, the magnetic resonance peak at 42 meV is 
centered at the AF wave vector and displays a Gaus- 
sian lineshape on top of a sloping background (Fig. [Hb). 
A similar signal can be observed at 40 meV, whereas 
at slightly lower energy, 38 meV, the magnetic response 
weakens and starts to broaden. Although the spin re- 
sponse is not clearly peaked at an incommensurate wave 
vector, the lineshape of the signal at this energy can be 
fitted to a double-peak structure. On decreasing the en- 
ergy transfer further, no sizeable magnetic response can 
be detected at 32 meV. Nonetheless, within the error bars 
one cannot exclude a weak magnetic response of magni- 
tude at most half of the one measured at 38 meV. It is 
worth noting that below 36 meV, the large nuclear back- 
ground makes the detection of a magnetic response (if 
any) particularly difficult. The broadening of the mag- 
netic response is not observed exclusively below the en- 
ergy of the magnetic resonance peak. At higher energy 
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FIG. 3: Dispersion of the resonant magnetic excitations de- 
duced from constant energy scans: a) along the (130) direc- 
tion (the full lines are described in the text), b) along the 
(110) direction. Horizontal and vertical error bars stand for 
the half width at half maximum of the signal and the energy 
resolution, respectively. 

(scans at 46 and 50 meV in Fig. [21b), the magnetic re- 
sponse also broadens, but in contrast to the measure- 
ments at 38 meV, the magnitude of the signal remains 
comparable with that of the magnetic resonance. While 
at 54 meV a magnetic response with a double peak pro- 
file can be observed, at slightly higher energy, 58 meV, 
the magnetic response seems to be weak and commensu- 
rate, suggesting that the spin excitations may move back 
to the AF wave vector. The broadening of the magnetic 
response above and below E r =A2 meV is consistent the 
differential energy scans performed away from the AF 
wave vector at Q=(0. 52, 1.56,4.7) and Q=(0.54,l. 62,4.7) 
(Fig. [TJb). Indeed, at these wave vectors, the magnetic 
scattering is widely spread in energy, but shows a slight 
dip at 42 meV. The combination of energy scans and 
constant-energy scans along the (130) direction is consis- 
tent with the existence of resonant spin excitations dis- 
persing upward and downward, as previously observed in 
the YBCO system. Guided by the data in YBCO±^ii, we 
then mo del both dispersions of F ig. [3]a by the relations: 
E± = \f{E r ) 2 ± (a±(q - <7af)) 2 where (-) and (+) cor- 
respond to the downward and upward dispersions, re- 
spectively. One obtains a_ = 130 meV.A and a + = 200 
meV.A in good agreement with the dispersions in YBCO 
obtained for similar doping level o 13 ! 14 ' 15 . 

Along the (110) direction, the evolution of the mo- 
mentum distribution of the magnetic signal as a function 
of the energy exhibits the same trends as those observed 
along the (130) direction, with a broadening at both high 
and low energies (46 meV and 40 meV in Fig. [51 a). A 
slight double-maximum structure also shows up in the en- 
ergy scan at Q=(0.425,0.425,14) (Fig.[TJa). But a closer 
inspection of the data reveals marked differences. First, 
the signal begins to broaden already at 40 meV, and be- 
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FIG. 4: a) Temperature dependencies at different wave vec- 
tors and energies, b) Differences between constant energy 
scans at 10 K and 100K. Scans are performed at 40 meV 
around the AF wave vector along 3 different directions: (100), 
(130), (110). The label corresponds to the angle between the 
scanning direction of the (100) direction (see Fig. [6jf ) . Data 
are normalized so that the intensities at the AF wave vector 
are the same. 



comes almost twice as broad at 38 meV. At this energy, 
the lineshape of the magnetic response is now rather dif- 
ferent, with a more complex structure that can be qual- 
itatively described in terms of a central peak at the AF 
wave vector surrounded by two satellites of similar mag- 
nitude at Q=(-0.5 ± 5, -0.5 ± S, 14) with 5=0.2. As a 
result of this expansion, a magnetic signal can still be 
observed in differential energy scans at Q=(0.35,0.35,14) 
(Fig. [Ua), i.e at a large distance in reciprocal space from 
the AF wave vector. 

As emphasized above, one of the hallmarks of the res- 
onant spin excitations in SC cuprates is their peculiar 
temperature dependence 1 -. In this respect, the behav- 
ior of Bi2212 is very similar to the phenomenology re- 
ported in YBCO. The magnetic resonance peak at 42 
meV disappears steeply at T c , as can be seen in Fig. |4ja 
at 42 meV. The temperature dependence measured at 
Q=(0.35,0.35,14) and 40 meV suggests a similar change 
at T c , indicating that the differential signal at this wave 
vector is also of magnetic resonant type. At the back- 
ground position, Q=(0.25,0.25,14), the T-dependence 
does not show any indication of a magnetic signal. 

We now discuss some general aspects of the spin excita- 
tions of Bi2212. While the overall layout of the magnetic 
spectrum of Bi2212 is quite similar to that of YBCO, 
well-defined incommensurate magnetic peaks were not 
observed in Bi2212. Several reasons can be put forward 
to explain this difference. Obviously, the weak signal- 
to-noise ratio of the measurements in Bi2212 limits the 
possibility to see details of the q-dependence. In partic- 
ular, we had to work with a broad q-resolution in order 
to pick up enough intensity. (Attempts with improved 
q-resolution were not successful because of insufficient in- 
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FIG. 5: Differences between constant energy scans at 46 meV 
performed between : a) 10 K and 300K, b) 100 K and 300K. 
Measurements were carried out on the spectrometer IN8. The 
solid line corresponds to the fit of the data by a double Gaus- 
sian functional form, on top of a sloping background (dashed 
line). The background at the AF wave vector is set to to 
get rid of the the variation of the thermal enhancement of the 
nuclear background and to allow a direct comparison of both 
results. 



tensity.) However, scans with similar resolution and scat- 
tering plane 13 did reveal incommensurate peaks below E r 
in YBCO. Therefore, the origin of the difference between 
both systems is not purely instrumental, but at least in 
part intrinsic to Bi2212. It was previously noticed 1 - that 
energy scans in Bi2212 exhibit a width larger than the 
energy resolution for all doping levels, as shown in Table 
HI whereas the energy width of the odd resonance peak in 
YBCO is resolution-limited at optimal dopin g 14 ' 15 . The 
intrinsic width of the mode blurs the details of the mode 
dispersion in Bi2212. If we nevertheless try to fit the 
data away from E r with a double peak structure, one ob- 
tains dispersive excitations both below and above E r as 
shown in Fig. [3l as previously reported in YBCO 13 ' 14 ' 15 . 
The overall energy and momentum dependences of reso- 
nant spin excitations then exhibit the typical hour-glass 
lineshape characterizing the resonant mode dispersion in 
YBCO, as evidenced by the color maps of Figs. [S]d- 
e deduced from our data along the two q-directions of 
Fig. [2 At E r —42 meV, the magnetic signal shrinks 
around qAF- Further, below E r , additional magnetic ex- 
citations appear to develop in a momentum region far 
from qAF- As shown by color maps of Figs. [6]d-e, the 
spectral weight of these excitations is maximum along 
the diagonal direction, leading to an anisotropy that we 
discuss in section [V] 

As pointed out at the beginning of this section, we 
focus on the enhancement of the magnetic response in 
the SC state. The main reason is that previous INS 
measurements in optimally doped Bi2212 10 show that 
normal state excitations are not measurable around 40 
meV. These early measurements are in agreement with 
the results in YBCO for the same doping level 5 -^. In 
YBCO^, a systematic study as a function of hole doping 
indicates that the normal-state AF fluctuations weaken 
continuously with increasing doping level and fall below 
the detection limit above optimal doping^. Note that 
the fact that the magnetic fluctuations are not sizeable 
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above T c does not necessarily mean that they are absent, 
but it rather suggests that they are much weaker and/or 
broader in the normal state than in the SC state. For 
instance in weakly overdoped YBCO, the typical mag- 
nitude of the spin fluctuations left in the normal state 
at the resonance energy is estimated to be one order of 
magnitude weaker that magnetic resonance peak in the 
SC stated. In our slightly overdoped Bi2212 sample, the 
upper limit on the magnitude of the magnetic signal left 
in the normal state extracted from our data is about 1 /4 
of the signal in the SC state at 46 meV (Fig. [5]). This 
is consistent with a previous report at 43 meV in an op- 
timally doped Bi2212 sample^. 



IV. FINGERPRINTS OF THE 
ELECTRON-HOLE SPIN FLIP CONTINUUM 

In order to elucidate the relationship between the res- 
onant spin excitations in the SC state and the gapped 
Stoner continuum, we have computed the threshold of 
the continuum in the odd channel for the two main di- 
rections along which most of the INS measurements were 
carried out (see Fig. [Slb-c): (130) and (110). One has 
to pay attention to the bilayer structure that affects 
both spin and charge properties. The motion of elec- 
tron between the two layers in a bilayer unit leads to the 
formation of anti-bonding (a) and bonding (b) states, 
and consequently to a splitting of the Fermi surface, as 
shown in Fig. [5Ja. The odd neutron scattering chan- 
nel originates from electronic interband spin flip excita- 
tions, and the threshold of the Stoner continuum in the 
SC state is defined as the minimum [E% + E%. + ], where 

E^ b — \J £j?' 62 + A 2 , is the electronic dispersion relation 
in the SC state. Afc = A m (cosfc K — cosfc y )/2 stands 
for the d-wave SC gap and the bare electronic dispersion 
relation is described using a simplified tight-binding ex- 
pression: = 2t(cosk x + cos k y ) — At' cos k x cos k y + 
2t" (cos 2k x — cos 2k y ) ± it^(cos k x — cos k y ) 2 — fi. For a 
nearly optimally doped Bi2212 sample with a T c of 87 K, 
we used the following parameters determined by ARPES: 
A m =35 meV 51 and {0.219,0.108,0.207,-0.95} for the 
parameters {£', t" , t±, /i} in units of t 41 . The chemical 
potential corresponds to a hole doping level determined 
from the SC transition temperature according to the phe- 
nomenological relationship of Tallon et a^. 

Before proceeding to the comparison of the locations 
of the gap in the Stoner continuum and the observed res- 
onant spin excitations, one needs to mention that there 
are two ways to determine the absolute scale of electronic 
parameters. The value of these parameters in absolute 
units (meV) strongly relies on the estimate of the value 
of the nearest-neighbor hopping parameter t, which in 
turn can be extracted from the Fermi velocity vf mea- 
sured in ARPES experiments. Along the nodal direction, 
vf — 2.0 eVA^S, is related to the bare Fermi velocity by 
vf = u f/(1 + ^) where A describes the electronic in- 



teractions related to the real part of the electronic self- 
energ y 29 ' 40 . In principle, v F can be determined by the 
band structure given by the Local Density Approxima- 
tion (LDA). In a self-consistent experimental procedure 
using the Kramers-Kronig transformation to determine 
the self-energy 4 ^, the energy scale t is deduced from an 
estimate of the bare Fermi velocity along the nodal di- 
rection, v° F ~ 4.0 eV.A 4 ^, yielding i=397 meV, in good 
agreement with LDA predictions. An alternative ap- 
proach to analyse the ARPES data essentially defines an 
effective band structure describing only low energy elec- 
tronic excitations (see e.gM-). The hopping parameter 
determined in this way is about a factor of two smaller: 
t ~ 200 meV. This does not significantly change the 
lineshape of the gapped Stoner continuum and can be 
ignored in this section. However, it is important for a 
quantitative description of the S=l mode origin (see sec- 
tion EE}. 

Figures [6]i-e show color maps of the magnetic inten- 
sity deduced from fits to the constant energy scans of 
Fig. [3 The threshold of the gapped Stoner continuum 
is superimposed on the magnetic signal. For the sake of 
simplicity, one can define two distinct areas below the 
continuum (Fig. [51b-c). In area I, the threshold of the 
continuum decreases away from the AF wave vector, ex- 
hibiting a dome-like shape. Area II corresponds to the 
remaining part of the gapped portion of the phase space. 
Along the (130) direction, the magnetic response in the 
SC state is confined to area I and appears to asymp- 
totically approach the threshold of the continuum. The 
signal vanishes below the detection limit in areas where 
continuum excitations are possible according to the cal- 
culation. 

Along the (110) direction, the situation is more com- 
plex. As shown in Fig. [5fc,e, area I is smaller than along 
(130), so that its spread in momentum becomes simi- 
lar to the q-width of the experimental resolution ellip- 
soid (Fig. [SJc). Moreover, the intrinsic energy width dis- 
cussed above scrambles the observed magnetic response, 
limiting an accurate determination of the exact number 
of branches of the magnetic dispersion and their pre- 
cise location. Some interesting observations can be made 
nonetheless. In particular, along (110) area II extends 
up to energies in the range 35-48 meV, higher than along 
(130). The extra magnetic signal along the (110) at 38 
meV may hence be attributable to secondary spin exci- 
tations in area II, in addition to those observed in area 
I. In slightly underdoped YBCOg.85, such a secondary 
magnetic contribution has been observed in area II, al- 
beit at higher energies around ~ 54 meV— . There it 
was also shown that the continuum gives rise to almost 
vertical silent bands 14 , where the intensity of collective 
magnetic modes is suddenly suppressed presumably due 
to the decay into elementary electron-hole spin flip ex- 
citations. This scenario was further confirmed theoret- 
ically within the spin exciton model^i. Interestingly in 
Bi2212 the scan along (110) at 38 meV (Fig. da) dis- 
plays a slight minima at the planar wave vectors q~(- 
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excitations directly measured by ARPES for the same 
doping level 4 ^, without adjustable parameters. This in- 
terpretation is also qualitatively consistent with the con- 
verse computation of the imaginary part of the dynamical 
susceptibility from ARPES data in the framework of the 
spin exciton scenario 45 ' 46 . These computations indeed 
indicate the existence of two magnetic resonant modes, 
located in area I and area II respectively. The detailed 
momentum shape of the magnetic spectra (Figs. [6]d,e) 
further suggests that the resonant spin excitation spec- 
tra exhibit a hourglass lineshape mainly confined to area 
I. Overall, these results are in good agreement with pre- 
vious studies carried out in YBCO syste m 14 i 18 i 19 i 21 , al- 
though the intrinsic energy width of the magnetic modes 
in Bi2212 obscures some of the features. 



ANISOTROPY OF THE MAGNETIC 
RESPONSE 





0.2 0.4 0.6 0.8 
Q=(H,H) 



FIG. 6: (color online) a) Fermi surfaces in nearly optimally 
doped Bi2212 deduced from ARPES measurement ;] 41 ' 5 1 . 
(red=bonding band, black=anti-bounding band). Threshold 
of the electron-hole spin flip continuum in the SC state along 
different directions: b) (130), c) (110). Below the threshold of 
the continuum, one distinguishes two distinct areas labelled 
respectively I and II as in^ (see text). The vertical dashed 
lines indicate the location of the energy scans which were per- 
formed. The projection of the resolution ellipsoid at 40 meV 
is also represented, d-e) Color maps showing the enhanced 
magnetic response in the SC state as deduced from the fit 
of the constant energy scan, in addition to the threshold of 
the continuum, f ) Area occupied by the continuum at 40 meV 
(black) . The red points indicate the extension of the magnetic 
as deduced from Fig.[4]b (see text) and the white crosses show 
where the temperature dependencies of the scattering inten- 
sity were measured. 



0.4,-0.4) and q~(-0.6,-0.6), which could be indicative of 
silent bands (see Fig. [Sic) in a location that coincides with 
the one expected based on the Fermi surface topology 41 . 

The comparison between the location of the observed 
magnetic excitation spectra and the momentum shape of 
the Stoner continuum of the d-wave superconductor thus 
shows that spin excitations are only present in the gapped 
regions of the continuum. We emphasize that the locus of 
the Stoner continuum was computed from the electronic 



Another interesting experimental feature is the in- 
plane anisotropy of the magnetic response. Figure [5]b 
shows scans performed at 40 meV around the AF wave 
vector along 3 different directions: (100), (130), and 
(110). The data are labelled as a function of the an- 
gle O between the direction of the scans and the (100) 
direction (see Fig. [6]f). Furthermore, the data are plot- 
ted as a function of the reduced distance to the AF wave 
vector (in units of A" 1 ), so that the q- widths of the sig- 
nals are directly comparable. In all directions, we fit 
the magnetic peak to a single Gaussian. From 8 = 0° 
where A g = 0.4 A" 1 (Full Width at Half Maximum) to 
= 18° where A q = 0.23 A -1 , the magnetic signal ex- 
hibits a Gaussian profile, with a reduction of its q-width 
at 6 = 18°. On the contrary at 6 = 45°, the q-width 
of the signal more than doubles, A q — 0.9 A -1 , indicat- 
ing a net anisotropy of the magnetic response along the 
(110) direction. Along that direction and at the reso- 
nance energy 42 meV, one obtains ~ 0.45 A" 1 for the in- 
trinsic q-width of the magnetic signal after deconvolution 
from the resolution function in agreement with previous 
report a 10 ' 11 . 

Concerning the origin of the observed anisotropy of 
the magnetic response along the diagonals, one needs to 
keep in mind that the structure of Bi2212 is not simply 
tetragonal. Indeed, Bi2212 is actually an orthorhombic 
system with an incommensurate distortion likely of com- 
posite type 5 ^, because the lattice parameters of the Cu02 
planes and those of the Bi02 planes do not match. The 
orthorhombic axes are along the diagonals of the square 
lattice: the sample growth direction, i.e. a* ortho , was 
kept perpendicular to the scattering plane. The direc- 
tion (110) in the tetragonal notation we have adopted 
in this article correspond to b* ortho . The incommensu- 
rate modulation is given by Q.21b* ortho + c* . Because of 
this modulation, the Bi2212 samples are not twinned, al- 
though the in-plane lattice parameters are not very dif- 
ferent. It is therefore conceivable that the observed mag- 
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netic anisotropy along b* ortho is related to this ID struc- 
tural anisotropy. It is interesting to remark here that in 
the YBCO system, the magnetic response also exhibits 
a ID-like anisotropy of the magnetic intensity, which is 
maximal along a particular in-plane direction , the (100) 
direction perpendicular to the CuO chains^. As the CuO 
chains in YBCO, the Bi02 planes in the Bi2212 system 
play the role of charge reservoir. We therefore cannot 
exclude that the diagonal anisotropy is a feedback of the 
structural distortion on the magnetic properties. How- 
ever, there are differences in the magnetic anisotropy of 
both systems. In the YBCO system, the anisotropy near 
optimal doping is mostly related to the intensity 2 ^.. Here, 
the q-extension of the magnetic signal appears to be dif- 
ferent along the two directions. Further INS experiments 
in which the two in-plane directions a* ortho and b* ortho are 
studied under the same resolution conditions, following 
prior work in YBC02&, will be needed to shed light on 
the origin of the anisotropy of the magnetic response in 
Bi2212. 

We report in Fig. Of the momentum dependence of the 
threshold of the continuum at 40 meV. The area I cen- 
tered at the AF wave vector exhibits a diamond shape, 
whereas the area II splits into 8 segments (4 along the di- 
agonals and 2 along a* and b* ) . The red points in the fig- 
ure indicate the total momentum expansion of the mag- 
netic responses (roughly twice the full width of the signal 
at half maximum), as deduced from Fig. [4]b. The data 
have been symmetrized to account for the fourfold sym- 
metry of the Cu02 plane although, as we have remarked 
above, it is not known whether the magnetic signal re- 
spects the square lattice symmetry. Anyway, one notices 
in Fig. Of that the magnetic signal matches the area I 
along (100) but extends into area II along the diagonals. 
However, the intensity of the magnetic signal in area II 
along the (110) direction is of the same order as the er- 
ror bars in the measurement performed along the (100) 
direction. We therefore cannot rule out the existence of 
a weak response in area II also along the (100) direction, 
but this extra magnetic response (if any) vanishes along 
the (130) direction where area II is absent. 

Finally, whatever the origin of the anisotropy, its ob- 
servation allows us to solve an old puzzling issue. In the 
initial work of Fong et al. in optimally doped Bi2213i£, it 
was found that the energy-integrated spectral weight at 
qAF of the magnetic resonance peak was similar in opti- 
mally doped Bi2212 and YBCO (~ 2(i%/f.u), whereas 
the local spin susceptibility (integrated in momentum 
space) was actually four times larger in Bi2212 than 
YBCO. This discrepancy came from the estimation of 
local spin susceptibility based on the assumption that 
the momentum distribution of the magnetic response was 
isotropic in Bi2212, as it is the case at qAF in YBCO. 
Since the measurement of the momentum width of the 
magnetic signal along the (110) direction of Bi2212 was 
found to be twice as large as in YBCO, the local suscepti- 
bility was estimated to be four times larger. Considering 
Fig. Of, the momentum shape at the resonance energy 



found in optimally doped YBCO typically corresponds 
to only area I, and extra scattering in area II occurs at 
larger energy^. In optimally doped Bi2212, this extra 
magnetic scattering in area II is found around the same 
energy as the resonance peak, yielding a much broader 
peak in q-space. Therefore, the present study shows that 
the assumption of an isotropic magnetic response was 
not correct, leading to an overestimation of the local spin 
susceptibility. This implies that the q-integrated spectral 
weight is likely quite similar in both systems. 



VI. ORIGIN OF THE S = 1 COLLECTIVE 
MODE: RPA DESCRIPTION 

Based on the energy and momentum distribution of the 
resonant spin excitations in the SC state, we have shown 
above that our INS data exhibit fingerprints of the d- 
wave gapped Stoner continuum. This puts constraints 
on the origin of the 5 = 1 collective mode. However, 
as far as the electron-hole spin flip continuum is taken 
into account, both itinerant and localized spin models 
can, in principle, describe the data. Since the charge 
excitation spectrum and the SC gap are well-known in 
the Bi2212 family from a considerable amount of ARPES 
data, this system is the right candidate to test these sce- 
narios quantitatively. Here, we apply the spin-exciton 
mode l 29 i 30 ' 31 i 32 ' 33 i 34 using the measured band structure 
parameters and SC gap. Furthermore, within that frame- 
work, one can relate the energy width of the resonance 
peak to the spatial distribution of the SC gap as reported 
by Scanning Tunnelling Microscopy (STM)^iii. We then 
investigate the magnetic excitation dispersion along the 
two directions studied by INS in order to put this com- 
parison on a quantitative footing. 

In the spin-exciton scenario, the spin susceptibility 
takes an RPA-like forrrtS 9 - : 



x(q,w) 



X°(q,w) 



l-V q X°(q,w) 



(3) 



X° (q, w) stands for the standard non- interacting spin 
susceptibility of a superconducto r 29 i 30 i 31 i 32 i 33 i 34 . The in- 
teraction V^j that enhances the magnetic response can 
be interpreted as: (i) the on-site Coulomb repulsion on 
copper U in weak coupling models (U < 8t, 8t being 
of the order of the band width), (ii) the AF superex- 
change coupling — J(q) = —2J(cosq x +cosq y ) in strongly 
correlated models (U > 8t, such as in the t-t'-J model, 
with J = At 2 /U ~ 120 meV) 3 ^' 3 ^, (hi) an effective spin- 
fermion coupling g(q) in phenomenological model o 30 i 34 . 

As before, we limit our calculations to the odd excita- 
tions. In this channel x°(q, u>) involves two-particle exci- 
tations between bonding (b) and anti-bonding (a) states: 
X° = (Xab + X° ba )/ 2 - To calculate x°(q,w), we use a 
SC gap with d x 2 +y i symmetry and a normal-state tight 
binding dispersion previously defined in section IIV1 In 
order to facilitate the numerical calculations, we used a 
small damping parameter of 2 meV, which is smaller than 
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FIG. 7: Resonant magnetic modes in nearly optimally doped 
Bi 2 Sr 2 CaCu208+a (T c =87 K) measured at Q = (0.5,0.5, 14) 
(as in Fig. [T]a). The solid lines correspond to computed 
imaginary part of the susceptibility convoluted with the en- 
ergy resolution function: the green line was obtained using a 
a single d-wave SC gap (RPA 1) and the red one with a distri- 
bution of d-wave SC gaps (RPA 2). In the latter case, the SC 
gap average is A m =35 meV with a full width at maximum of 
14 meV. 



the energy width of the INS resolution. Finally, we use 
the following phcnomcnological form for the interaction 
= U e jf — 2J e ff(cosq x + cosg y ) 34 . This allows us to 
capture the different physical origins of the interaction. 

For any given wave vector q, a sharp magnetic ex- 
citation shows up at an energy £l r when the conditions 
1 - F q i?ex°(q, fi r ) = and Imx°(q, fl r ) = are fulfilled. 
For a given magnitude of the SC gap amplitude A m , the 
measured energy position of the odd magnetic resonance 
peak at qAF = (7r, t) constrains the sum of the interaction 
parameters Vq AF — U e ff + 4J e //; whereas the shape of 
the dispersion is controlled by their ratio J e ///i7 e //. 

Let us start by studying the odd magnetic resonance 
peak at the planar AF wave vector qAF- In Fig- 13 we 
compare the measured resonant peak energy scan at qAF 
and the spin excitation spectrum computed with Eq. [3] 
and convoluted with the energy resolution of the spec- 
trometer. For a slightly overdoped sample with T C =87K, 
A TO is set to 35 meV, yielding a threshold of the Stoner 
continuum of cj c =62 meV at qAF- To obtain a reso- 
nance energy at the proper energy, f2 r =42 meV, well be- 
low lo c , one needs to adjust the interaction Vq AF =l§7§ 
meV, a value significantly exceeding the hopping param- 
eter i=397 meV. A systematic investigation yields a ratio 
Vq AF /t = 2.7, independent of t in the range 200-400 meV, 
in agreement with Ref. [3J. 

At qAF, the computed spin excitation spectrum is 
dominated by the spin exciton, whose contribution to 
the imaginary part of the dynamical spin susceptibility 
reads: Imx(q_AF,^) — W r 5(u> — fi r ) for uj > 0. The 
measured magnetic resonance peak should be given by 



the convolution product of this 5-function by the Gaus- 
sian resolution function of the instrument: it shquld take 
the form, /mx( c UF,w) = W r exp(— 4 In 2 ^ 4^ ) where 

<7 W ~ 6 meV is the energy resolution (as shown by the 
green line (RPA 1) in Fig. [7]). Clearly, the energy width 
of the resonant mode is significantly broader than the 
expected spectrum, showing that there is an additional 
contribution to the broadening. This discrepancy is not 
specific to our sample. Table fl] shows the corresponding 
energy position and energy width of all Bi2212 samples 
investigated so far by neutron scattering. For all doping 
levels, the measured resonant mode is broader than the 
calculated resolution-limited RPA peak. For a long time 
now, it has been argued 1 - that this broadening might be 
related to the electronic inhomogeneity of the Bi2212 sys- 
tem, which also manifests itself in the SC gap distribution 
observed in the real space by ST M 42 i 43 i 44 . Interestingly, 
when going into the overdoped regime, the broadening 
of the INS magnetic resonance peak 4 ^ and the SC gap 
distribution 4 ^ are reduced simultaneously (see Tablc|T]for 
the resonance peak). 

Based on this observation, we have tested the effect of a 
spatial distribution of SC gap on the spin excitation spec- 
trum within the spin exciton scenario. Spatial variations 
of the SC gap at the surface of Bi2212 were reported by 
several group a 42 ' 43 i 44 . The spatially averaged value of the 
superconducting gap, A m , measured by STM is consis- 
tent with other spectroscopic techniques such as ARPES. 
The FWHM of the gap distribution, a a , decreases with 
increasing doping level (Fig. [9]), from around 15 meV 
at optimum doping to nearly 7 meV for an overdoped 
sample with a gap average of A m =23 meV— . In our cal- 
culation, we use the gap distribution found by STM at 
optimal doping, and we further assume that the inter- 
action Vq is independent of the local SC gap. Since the 
local gap amplitude varies, the resonance peak position 
is different in different patches of the sample. We com- 
pute the susceptibility x(q, u>) in each patch character- 
ized by a different SC gap value, labelled now x(q, uj, A). 
x(q, ijJ, A) is computed for seventeen different SC gaps 
with A from 20 to 50 meV. The gap distribution is im- 



10 



30 



25 



> 

CD 

£Z 
O 

3 

CO 
T3 

Q. 
CO 

CD 



20 



15 



10 




0.1 



0.15 
P 



0.2 



0.25 



FIG. 9: Width (FWHM) of the SC gap distribution measured 
by ST M 42 ' 43 ' 44 and width of the magnetic resonance peak in 
INS dataia^. 



plemented using the histogram shown in Fig. [51 which 
approximates a Gaussian distribution of the SC gap with 
a width a a- The full susceptibility, x(q, u), is given by 
the sum over all local susceptibilities x(q, w , A): 



x(q,w) 



x(q, u>, A) exp(— 4 In 2 



(A - A m f 



)dA (4) 



The total spin susceptibility is further convoluted by the 
Gaussian resolution function in order to fit the data. 

Using this procedure, we can simulate the broadening 
of the magnetic resonance peak by adjusting the width 
of the Gaussian gap distribution, a a (Eq. Figure [7| 
(fit: RPA 2) shows that the enhancement of the mag- 
netic response in the SC state at the AF wave vector can 
be well accounted for by an intrinsic SC gap distribution 
with cta = 14 meV. For an overdoped sample with T c =70 
K4 8 -, the same analysis gives cta=5 meV. It is striking 
that the gap distribution deduced from the fit of the 
magnetic resonance peak matches the one reported by 
STM for both doping levels (Fig. [9]). In the framework 
of spin exciton scenario, one can thus give a clear ex- 
planation of the energy width of the magnetic resonance 
excitation. Note that such energy broadening is absent 
in the YBCO family compounds. This shows that the 
electronic inhomogeneity is not generic to the cuprates. 
Similar conclusions about a better homogeneity in YBCO 
were reached by measuring the quasi-particle lifetimes 56 
or the 89 Y nuclear magnetic resonance linewidths 57 . 

Next, we address the consequences of the SC gap dis- 
tribution for the q-dependence of magnetic excitations, 
still within the spin- exciton model. We performed the 
calculation in the q — lu range covered by our experi- 
ments, i.e. along both the (110) and (130) directions 
and between 30 and 60 meV. Using the same set of pa- 
rameters, we additionally adjust the ratio J e ///£/ e // to 



obtain the best agreement between the experimental q- 
dependence (Figs. [(J]d,e) and the computed maps of Figs. 
rTTJl This corresponds to the case when J e ff/U e ff << 1. 
We compute the maps of FigfTTJlwith J e ff /U e ff = 0.025 
for both Q directions. As for the energy width at the 
AF wave vector, the model does not describe the results 
in the absence of a SC gap distribution (Fig. [TUJa and 
d). One obtains a rather sharp mode dispersion down- 
ward as observed in YBCOi^. In contrast (Fig. |T0]b and 
c), the addition of a SC gap distribution (correspond- 
ing to the STM data) allows us to describe the main 
features of the measured magnetic excitation spectrum 
(Fig. [lOjc and f). For the (130) direction, the experi- 
mental and computed mappings are quite similar. For 
the (110) direction, both mappings are consistent, even 
if below 42 meV the computed spectrum fails to describe 
the q-broadening of the signal. We found that it is not 
possible to obtain a signal around 38-40 meV in area II, 
whatever the q-dependence of the interaction used. 

To summarize, using the band structure and the SC 
gap known from ARPES measurement s 41 ^ 1 , the calcula- 
tion of the magnetic excitations in an itinerant spin ap- 
proach can reproduce most of the characteristic features 
of the INS spectrum measured the spin excitation spec- 
trum in the SC state of Bi2212. Further, the observed 
energy width of the resonant spin excitations in Bi2212 
is naturally explained within this model as a result of 
the gap distribution reported by STM^^di. The mea- 
sured q-dependence of the magnetic excitations can be 
also captured by a broadening of dispersive excitations 
due to the SC gap distribution. 

Finally, we comment on a few points about the nature 
of the interaction, that should provide an indication on 
the Hamiltonian needed to obtain a spin exciton. 

First, one finds a ratio Vq AF /t = 2.7 whatever the 
used band structure. We consider the energy range for 
t which is typically 200-400 meV. Two opposite limits 
can then be discussed: (i) a weak coupling approach, 
where the band structure is not far from the LDA one, 
i.e. t ~ 400 meV, then the interaction is an effective on- 
site Coulomb repulsion, (ii) a strong coupling approach, 
where one chooses an effective band width, i.e. t ~ 200 
meV, reduced from the LDA value due to the strong elec- 
tronic correlations and with an interaction given by the 
AF superexchange interaction J (in principle measured 
in the AF insulating state). In the first limit, using £=395 
meV, the effective interaction U qAF = 1070 meV can be 
directly compared to the band width, W = 8t. One finds 
Vq AF /W « 0.35, which clearly belongs to the weak cou- 
pling approach 6 ^ with an effective on-site Coulomb inter- 
action on copper. In the other limit, t is around 200 meV 
and then V^ F =540 meV. This value can then be reduced 
to ~ 4 J as expected in the t — t' — J model 32 , yielding 
J = 135 meV. In the insulating state, J is usually re- 
lated to the on-site Coulomb interaction as J = 4t 2 /U, 
yielding [7=1185 meV. Surprisingly, the ratio between 
the interaction and the band width is found to be quite 
small as U/W ~ 0.74, i.e. in an intermediate regime not 
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FIG. 10: Mappings of the imaginary part of the dynamical 
spin susceptibility /mx(q,i<j): ( a_c ) along the (130) direction, 
(d-f) along the (110) direction. /mx(q,w) was computed us- 
ing a spin itinerant model (see text) assuming: (a,d) a unique 
SC gap for the entire system or (b,e) the spatial distribution 
of SC gaps reported in Fig. [5] In order to enable a immediate 
comparison of the computed mappings of Jmx(q,u) and the 
INS measured ones, the experimental mappings of Fig. [6jd- 
e are also reported (c,f). In absence of a calibration of the 
INS data in absolute units, mappings are scaled so that the 
maximum intensity corresponds to 100 counts. 



far from the weak coupling side. 

Next, we found that is rather weakly dependent on 
the wave vector, since J e ff/U e ff « 1. It is worth to 
note that in case of a momentum independent interac- 
tion, i.e. J e ff = 0, Eq. [3] is unstable towards a spin- 
density wave state using the band structure observed by 
ARPES in Bi2212. In the other limit, still taking the 
same electronic parameters but only a q-dependent in- 
teraction, i.e. U e ff = 0, one observes a resonance peak 
which first disperses upward and then downward when 
approaching the electronic continuum delimiting the area 
I and II in Fig. M (i.e. the silent bands). Such an "M"- 
shaped dispersion is not observed in Bi2212, although 
the SC gap distribution renders the experimental situa- 
tion too complex to rule this out entirely. Further, it is 
shown in Ref. [3l] that the exact magnetic mode disper- 
sion is extremely sensitive to the detailed tight-binding 
parametrization of the band structure. Therefore, de- 



pending on details of the band structure, one can easily 
move from a weakly to a strongly q-dependent interac- 
tion. 

Lastly, an important issue with the spin exciton model 
is related to the apparent observed "X" -shaped or hour- 
glass dispersion. Within the RPA approach it is difficult 
to reproduce downward- and upward-dispersing branches 
merging around f2 r . This is because in the simple RPA 
susceptibility of Eq. G2 there is only a single pole for a 
given wave vector below the electronic continuum, espe- 
cially in area I. The high-energy part can be understood 
by a pole condition of the RPA susceptibility below the 
electronic continuum, but only in area I I 14 ' 34 , without 
direct continuity with the mode at (7r,7r) and fl r . This 
limitation of the spin-exciton model may indicate strong- 
coupling effects not capture by the RPA formalism. The 



capability of alternative approximation scheme; 
reproduce this feature remains to be elucidated. 



,58.5!) 
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VII. CONCLUDING REMARKS 

We have shown that the spin excitation spectrum in 
the odd channel of Bi2212 is consistent with the hour 
glass or "X" -shaped dispersion previously reported in 
YBCO for a similar doping leve l 14 ' 15 . Most aspects of 
the momentum and energy dependence can be ascribed 
to the existence of a S = 1 collective mode, a spin- 
exciton. The study further reveals that the imaginary 
part of the dynamical magnetic susceptibility is strongly 
enhanced below the threshold of the electron-hole spin 
flip continuum that can be derived from ARPES mea- 
surements performed on the same system. This obser- 
vation suggests that the S = 1 collective mode decays 
into elementary electron-hole spin flip excitations when 
it enters the continuum. Likewise, the influence of the 
gapped Stoner continuum is confirmed by a recent study 
of the doping dependence of the characteristic energies 
and spectral weights of the resonance peaks in both odd 
and even channels^. In agreement with a similar study 
carried out over a wide doping range in the (Y,Ca)BCO 
system 5 - 5 -, the spectral weight of both AF resonance peaks 
is proportional to their reduced binding energy with re- 
spect to the continuum. In our Bi2212 sample (T c =87 
K), the estimate of the threshold of the continuum at the 
AF wave vector from INS data, oj c = 63 meV— , agrees 
remarkably well with the value deduced from ARPES 
measurements. 

As mentioned in the introduction, localized-spin mod- 
els based on the formation of stripe arrays of spins and 
charges have been also developed to account for the spin 
dynamics in high-T c cuprates. The predictions of these 
models can also be compared to our INS data obtained 
below T c in nearly optimally doped Bi2212 sample. How- 
ever, in most of these models the effect of the super- 
conductivity has not been addressed theoreticall y 25 ' 26 ' 27 . 
One stripe model that explicitly considers the presence 
of superconductivity predicts only minor effects, in con- 



12 



trast to our observations 61 . Most of this theoretical work 
is based on spin-only models, that is, they consider the 
magnetic response of localized Cu spins, but typically ig- 
nore the charge degree of freedom. However, it has been 
argued^ that when the spin-gap is large and close to 
the saddle-point energy, much of the magnetic spectral 
weight redistributed below T c is accumulated at the sad- 
dle point, yielding the strong commensurate resonance 
peak. In principle, this argument should apply to the 
Bi2212 system where the spin-gap is around 32 meV (Fig 
[T]). However, our data provide evidence that the key fac- 
tor to understand the spin dynamics in the SC state is 
the momentum and energy shape of the d-wave electron- 
hole Stoner continuum. In this picture, the spin-gap is 
ascribed to the energy where the collective mode merges 
into the continuum. Even in a stripe scenario, a gapped 
continuum delimiting the spin excitations must also be 
present. Our study tells us that this continuum has to 
be similar to the one obtained starting from uniform 
2D Fermi liquid theories. One more general grounds, 
dual (itinerant/localized) spin models incorporating the 
gapped Stoner continuu m 35 ! 36 ' 37 ! 38 may also be consis- 
tent with our data. 

Two features of the magnetic dynamics of Bi2212 are 
at variance with the YBCO system. First, the resonance 
peak in Bi2212 consistently exhibits an intrinsic energy 
width 1 — ' 11 ' 48 . Within the spin-exciton model, we were 
able to relate this width to the SC gap distribution ob- 
served by STM^JZdl (Fig. [9]). Second, we observed 
in Bi2212 an anisotropy along the diagonal(s) (110) di- 
rection, whereas in YBCO the spin excitation spectrum 
exhibits an energy dependent ID-like anisotropy^ with 
maximum spectral weight along a* . Whatever the origin 
of these anisotropics, this difference between YBCO and 
Bi2212 system suggests that the specific form of the in- 
plane anisotropy of the spin excitations may depend on 



structural details of individual compounds. 

In summary, we have determined the detailed mo- 
mentum dependence of the resonant spin excitations in 
the SC state of a nearly optimally doped Bi2212 sam- 
ple (T c —87 K. The salient features of this spectrum are 
well described in an itinerant-electron approach. Inter- 
estingly, the intrinsic energy width and its doping de- 
pendence are naturally explained by considering the SC 
gap distribution as measured by STM. The global mo- 
mentum shape of the measured magnetic excitations is 
also correctly described within the spin-exciton model 
in the presence of the same SC gap distribution. The 
spin excitation spectrum in Bi2212 is delimited by the 
gapped Stoner continuum in the d-wave superconducting 
state, which can be directly inferred from ARPES data 
on the same system. Together with prior observations 
in the YBCO system 1 ^, this underscores the influence 
of the Stoner continuum as an important ingredient in 
the description of the spin dynamics in superconducting 
cuprates with large T c . 
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